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ABSTRACT 


The combined effects of chemical reaction, radiation, Dufour and Soret effects on Casson MHD fluid flow over a vertical plate 
with heat source/sink has been studied with constant velocity, variable temperature and concentration. The governing equations 
of the flow have been solved by using the finite difference method. The effects of various parameters on velocity, temperature 
and concentration fields are studied with the help of graphs. It is found that velocity decreased with the increase of Casson pa¬ 
rameter, Hartmann number, thermal conductivity, Schmidt number, Chemical parameter and decrease of Soret effect, Thermal 
Grashof number and Mass Grashof number. Concentration found to be decreasing with the decrease of Soret number and 
increasing chemical parameter. Temperature decreased when thermal radiation is increased and heat source is decreased. 
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INTRODUCTION 

Casson fluid exhibits yield stress. It is well known that Cas¬ 
son fluid is a shear thinning liquid which is assumed to have 
an infinite viscosity at zero rate of shear, a yield stress below 
which no flow occurs, and a zero viscosity at an infinite rate 
of shear, i.e., if a shear stress less than the yield stress is ap¬ 
plied to the fluid, it behaves like a solid, whereas if a shear 
stress greater than yield stress is applied, it starts to move. 
Casson fluids like Honey, human blood, jelly, sauces, etc., 
are having significant importance. Human blood can also be 
treated as Casson fluid. 

In recent days, Casson fluids have occupied a prominence 
role in various fields like biological, chemical, medical and 
engineering applications, in industrial environment, Various 
researchers have carried out their research work in this di¬ 
rection enlightening the effects of various parameters on 
the flow of the Casson fluid.Mustafa and Khanstudied the 
MHD flow and heat transfer of Casson nanofluid over a non- 
linearly stretching sheetwith non-linear temperature distribu¬ 


tion.The heat and mass transfer characteristics of a magneto 
hydrodynamic Casson fluid flow on a parallel plate channel 
having stretching walls, when subjected to a uniform trans¬ 
verse magnetic field are analyzed by Sarojamma et al. Swati 
et al numerically examined the boundary layer flow due to 
an exponentially stretching surface in the presence of an ap¬ 
plied magnetic field. Casson fluid model is used to charac¬ 
terize the non-Newtonian fluid behavior. The flow issubject- 
ed to suction/blowing at the surface. Analysis is carried out 
in presence of thermal radiation and prescribed surfaceheat 
flux. Prakash et al studied the MHD two-dimensional bound¬ 
ary layer flow of a Casson fluid in the presence of chemical 
reaction and thermal radiationaccompanied by heat and mass 
transfer towards an exponentially stretching sheet. 

Animasaun et al. studied the thermal conductivity and mo¬ 
tion of temperature dependent plastic dynamic viscosityof 
steady incompressible laminar free convective MHD Cas¬ 
son fluid flow over an exponentially stretching surface with 
suction and exponentially decaying internal heat generation. 
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The effects of thermal radiation, suction/blowing, viscous 
dissipation, heat source/sink and chemical reaction onbound- 
ary layer flow of a non-Newtonian Casson fluid in presence 
of heat and mass transfer towards a porous exponentially 
stretching sheet with velocity slip and thermal slipconditions 
were analyzed by Saidulu and Lakshmi. The steady flow and 
heat transfer of Casson fluid from a permeable horizontal 
cylinder in the presence of slip condition in a non-Darcy po¬ 
rous medium is analyzed by maintaining the cylinder surface 
at a constant temperatureis studied by Prasad et al. Shehzad 
et al derived a series solution for the effect of mass transfer 
in the MHD flow of a Casson fluid over a porous stretching 
sheet in the presence of a chemical reaction. Casson fluid 
flow over a vertical porous surface with chemical reaction 
in the presence of hydromagnetic field has been studied by 
Arthur et al. The heat and mass transfer effect in a boundary 
layer flow of an electrically conducting viscous fluid sub¬ 
ject to transverse magnetic field past a moving vertical plate 
through porous medium in the presence of heat source and 
chemical reaction is analysed by Tripathy et al. 

The aim of the present paper is to analyze the combined ef¬ 
fects of chemical, radiation, Dufour and Soret parameters 
on Casson MHD fluid flow over a vertical plate with heat 
source/ sink has been studied with constant velocity, variable 
temperature and concentration. The governing equations 
of the flow have been solved by using the finite difference 
method. The effects of varies parameter on velocity, tem¬ 
perature and concentration fields are studied with the help 
of graphs. 

Formulation of the problem 

Consider a two-dimensional steady Casson fluid of incom¬ 
pressible, viscous, electrically conducting fluid over a verti¬ 
cal plate moving with constant velocity with radiation and 
chemical reaction in the presence of Soret and Dufour is 
considered. Consider the flow of an incompressible viscous 
fluid passing a flat sheet coinciding with plane y = O.We se¬ 
lect the Cartesian coordinate system such that the x - axis 
be taken parallel to the surface and y is perpendicular to the 
surface. The fluid occupies a half space y>0. The flow is 
subjected to a constant applied magnetic field B 0 in the y-di- 
rection. The magnetic Reynolds number is considered to be 
very small so that the induced magnetic field is negligible in 
comparison to the applied magnetic field. The surface tem¬ 
perature of the plate oscillates with small amplitude about a 
nonuniform mean temperature. The fluid is assumed to have 
constant properties except for the influence of the density 
variations with temperature and concentration which are 
considered only in the body force term. The temperature of 
the plate oscillates with little amplitude about a non-uniform 
temperature. 

The rheological equation of state for an isotropic flow of a 
Casson fluid [000]can be expressed as: 
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In Equation (1), n=e..e. p where eUs the (i , j) th component 
of deformation rate. This means that ^is the product of the 
component of deformation rate with itself. Also, n c is a criti¬ 
cal value of this product based on the non-Newtonian model, 
Mb is the plastic dynamic viscosity of the non-Newtonian 
fluid and P y is the yield stress of the fluid. The equations 
governed the unstudied boundary layer flow of the Casson 
fluid is 
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Equations (2),(3) and (4) refers Momentum Equation, En¬ 
ergy Equation and Species Equation respectively where u is 
the velocity of the fluid, (3 is Casson parameter, Q 0 is the 
heat source/sink parameter, D is the molecular diffusivity, k 
is thermal conductivity, C is mass concentration, t is time, n 
is the kinematics viscosity, g is the gravitational constant, (3 
and (3* are the thermal expansions of fluid and concentration, 
p is density, c p is the specific heat capacity at constant pres¬ 
sure, Dm Coefficient of mass diffusivity, kt Thermal diffusion 
ratio, Tm Mean fluid temperature, Too Free stream tempera¬ 
ture of the surrounding fluid, Coo Free stream concentration, 
T Fluid temperature, C Fluid concentration y is distance, qr is 
the radiative flux, (3 is the magnetic field, kr is the chemical 
reaction rate constant. 

R.H.S. of equation (2), second term is thermal heat 
effect,third term is thermal concentration effect, fourth term 
is magnetic effect, and second term is thermal buoyancy ef¬ 
fect. R.H.S. of equation (3) second term is thermal radiation 
flux and third term is thermal radiation and fourth term is 
Dufour effect. R.H.S. of equation (4), second term is chemi¬ 
cal reaction and third term Soret effect. Under the above as¬ 
sumptions the physical variables are functions of y and t. 
The boundary conditions for the velocity, temperature and 
concentration fields are: 

u = 0, r = Tl, c* = c: for all y ,f < 0 

u - U, T* -T* + e(t*-T*^ e nt , C* - C* + s(c*-Cfje nt at >> = 0 
>0, r ->0, C* —> 0 as 

Introducing the dimensionless quantities with thermal radia¬ 
tion flux gradient expressed and we assume that the tempera¬ 
ture differences within the flow are sufficiently small so that 
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* 4 

T can be expressed as a linear function of T after using 
Taylor’s series to expand T * 4 about the free stream tempera¬ 
ture t* and neglecting higher-order terms. This results in the 
following approximation. 

= -4 a o' (r; 4 - r 4 ) and T 4 = 4 T*T’ - 3r„* 4 ^ 


Method of Solution 

Equations (10)-(12) are linear partial differential equations 
and are to be solved with the initial and boundary conditions 
(13). In fact the exact solution is not possible for this set 
of equations and hence we solve these equations by finite- 
difference method. The equivalent finite difference schemes 
of equations for (10)-(12) are as follows: 


The following dimensionless quantities are introduced 
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The thermal radiation flux gradient may be expressed as fol¬ 
lows 
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Here, the suffix ‘i’ refer to y and ‘j’ to time. The mesh system 
is divided by taking Ay = 0.1. From the initial condition in 
(13), we have the following equivalent: 


Considering the temperature difference by assump¬ 
tion within the flow are sufficiently small such that 
r* 4 may be expressed as a linear function of the 
temperature. This is attained by expanding inr' 
Taylor’s series about T* 4 and ignoring higher orders 
terms. 

T* 4 =4T* 3 T* -3T* 4 (9) 


u(i,0) = 0, #(z,0) = 0,C(z,0) = 0 foralli (17) 

The boundary conditions from (13) are expressed in finite- 
difference form as follows 

u(0,j) = l, 0(0,j) = l +£e n ^ At , (p(0,j) = l + £e n ( J ^ At forallj 

(0,J) = 0, 0(O,j) = O, <p(0,j)=0 forallj ( 18 ) 


Substituting the dimensionless variables (7) into (2) 
to (4) and using equations (8) and (9), reduce to the 
following dimensionless form. 
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The corresponding boundary conditions are 


(Here i was taken as 200) 

The velocity at the end of time step viz, u(i, j+l)(i=l,200) 
is computed from (14) in terms of velocity, temperature and 
concentration at points on the earlier time-step. After that 0 
(i, j +1) is computed from (15) and then C (i, j +1) is com¬ 
puted from (16). The procedure is repeated until t = 0.5 (i.e. j 
= 500). During computation At was chosen as 0.001. 


Skin-friction: 

The skin-friction in non-dimensional form is given by the 
relation 


T 



, where T =-t- 

(V o 


u = 1, 0=l + ee nt , (p= l + £e nt at y = 0 
u = 1, 0=l + £e nt , (p= l + £e nt at y = 0 
u — ^ 0, 0 —> 0, (p — ^ 0 as y — > °° (13) 

Where G r is thermal Grashof number, P r is the prandtl num¬ 
ber, k is the chemical reaction parameter, R is the thermal 
radiation conduction number, M is Hartmann number, G is 
the mass Grashof number, Q is the heat source/sink param¬ 
eter and S is the Soret number. 

r 


Rate of heat transfer: 

The dimensionless rate of heat transfer in terms of Nusselt 
number is given by 


Nu = - 


\w_\ 

v d v J ,= 0 
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Rate of mass transfer: 

The dimensionless rate of mass transfer in terms of Sher¬ 
wood number is given by 


Sh = - 


r dC] 

V dy )y =o 


RESULTS AND DISCUSSION 

The effects of various parameters such as Casson 
Parameter “/?”, Thermal Grashof Number “G r ,„ 
Mass Grashof Number “G ”, Thermal Conductivity 
“K”, chemical reaction parameter k, Prandtl Number 
“P ” Schmidt Number “S ” Thermal radiation conduc- 

r c 

tion number “R”, Hartmann number M, Heat source/ 
sink parameter Q and Soret number S r on the veloc¬ 
ity, temperature and concentration fields are stud¬ 
ied numerically and represented the results through 
graphs. The influence of Casson parameter, Hart¬ 
mann number, thermal conductivity, Schmidt number 
and chemical reaction parameter on velocity is shown 
in the figures 1,5,6,8,9. From these figures, it is recog¬ 
nized that the velocity decreases with the increasing of 
these parameters respectively. The influence of Soret 
number, thermal Grashof number, the mass Grashof 
number on velocity are shown in the figures 2,3,4. It 
is clear that the velocity decreases with the decrease 
of these parameters respectively. The effect of Prandtl 
number on the velocity is shown in the figure 7. It 
is observed that the velocity decreases near the plate 
and increases far away the plate with the falling of 
the Prandtl number. The effect of Soret number on the 
concentration field is illustrated in figure 10. As the 
Soret number decreases the concentration is found to 
be decreasing. The effect of chemical reaction param¬ 
eter on the concentration field is illustrated in figure 
12. It is clear that the concentration is decreasing for 
increasing chemical reaction parameter. The effect of 
thermal radiation conduction number on concentra¬ 
tion field is shown in figure 11. It is noted that the 
concentration decreases near the plate and increases 
far away the plate with falling thermal radiation con¬ 
duction number. Figure 13 shows the variations of 
thermal radiation conduction number on temperature. 
It is found that the temperature is decreased when R 
is increased. Figure 14 shows the variation of the heat 
source sink parameter on temperature. It is clear that 
the temperature decreases when Q decreases. In Fig¬ 
ure 15, the skin friction decreases when thermal Gra¬ 
shof number increases and figure 16 shows that the 
skin friction decreases when Soret number increases. 



Figure 1: Effects of Casson Parameter/? on velocity. 



Figure 2: Effects of Soret Numbers,, on velocity. 



Figure 3: Effect of Thermal Grashof Number G r on velocity. 
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Figure 4: Effect of Mass Grashof Number G c on velocity. Fi 9 ure 7: Effect of Prandtl Number P r on Velocit V- 




Figure 5: Effect of Hartman NumberM on velocity. 



Figure 6: Effect of Thermal Conductivity K on Velocity. 


Figure 8: Effect of Schmidt Number S c on Velocity. 
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Figure 10: Effect of Soret Number P r on Velocity. 



Figure 11: Effect of Prandtl Number P r on Concentration. 



Y 

Figure 12: Effect of Chemical Reaction parameter^ on Con¬ 
centration. 



Y 

Figure 13: Effect of Thermal Radiation Conduction Number 
Ron Temperature. 



Y 

Figure 14: Effect of Heat Source /Sink Parameter Qon Tem¬ 
perature. 



Figure 15: Effect of Thermal Grashof Number G r on Skin 
Friction. 



Figure 16: Effect of Soret Number S r on Skin Friction. 
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CONCLUSION 

In this paper, from the study of the combined effects of chem¬ 
ical reaction, radiation, Dufour and Soret effects on Casson 
MHD fluid flow over a vertical plate with heat source / sink, 
it is concluded that 

• Velocity decreased with the increase of Casson param¬ 
eter, Hartmann number, thermal conductivity, Schmidt 
number, Chemical parameter and decrease of Soret 
effect, Thermal Grashof number and Mass Grashof 
number. 

• Concentration is found to be decreasing with the de¬ 
crease of Soret number and increasing chemical pa¬ 
rameter. 

• Temperature is decreased when thermal radiation is 
increased and heat source is decreased 

• Skin friction decreased when thermal Grashof number 
and Soret number are increased 
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